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There is ample evidence from rotation curves that dark matter halos around disk galaxies 
have nontrivial dynamics. Of particular signihcance are: a) the cored dark matter prohle 
of disk galaxies, b) correlations of the shape of rotation curves with baryonic properties, 
and c) Tully-Fisher relations. Dark matter halos around disk galaxies may have nontrivial 
dynamics if dark matter is strongly self interacting and dissipative. Multicomponent 
hidden sector dark matter featuring a massless ‘dark photon’ (from an unbroken dark 
U{1) gauge interaction) which kinetically mixes with the ordinary photon provides a 
concrete example of such dark matter. The kinetic mixing interaction facilitates halo 
heating by enabling ordinary supernovae to be a source of these ‘dark photons’. Dark 
matter halos can expand and contract in response to the heating and cooling processes, 
but for a sufficiently isolated halo could have evolved to a steady state or ‘equilibrium’ 
conhguration where heating and cooling rates locally balance. This dynamics allows the 
dark matter density prohle to be related to the distribution of ordinary supernovae in the 
disk of a given galaxy. In a previous paper a simple and predictive formula was derived 
encoding this relation. Here we improve on previous work by modelling the supernovae 
distribution via the measured UV and Ha huxes, and compare the resulting dark matter 
halo prohles with the rotation curve data for each dwarf galaxy in the LITTLE THINGS 
sample. The dissipative dark matter concept is further developed and some conclusions 
drawn. 
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1 Introduction 


Large scale structure (LSS) and measurements of the cosmic microwave background 
(CMB) provide strong evidence for the existence of nonbaryonic dark matter in the Uni¬ 
verse, e.g. [1, 2, 3, 4, 5, 6, 7, 8, 9]. Dark matter, whether collisionless or collisional, 
provides a compelling explanation for the measured anisotropies of the CMB and the 
observed large scale structure. On much smaller scales, measurements of the rotation 
curves of disk galaxies have also yielded very convincing evidence for the existence of 
nonbaryonic dark matter. Indeed, galaxy rotation curves are (typically) observed to be 
flat at the observed edge of the galaxy, in sharp contrast to the expected Keplerian decline 
[10, 11, 12, 13] (see also [14] and references therein). 

This asymptotic flatness is certainly an intriguing feature, yet rotation curves provide 
much more information. Detailed studies, especially those involving low surface brightness 
galaxies, e.g. [15, 16, 17] and gas rich dwarfs, e.g. [18, 19, 20], have found that the dark 
matter prohle is generally cored (see also the review [21]) and that the shape of rotation 
curves are correlated with baryonic properties. There are also local such correlations 
observed in specihc galaxies, such as the rotation curve ‘wiggle’ in NGC1560 [22, 23]. 
Similar conclusions have also been reached in high surface brightness spiral galaxies [24, 25, 
26]. See also [27, 28, 29, 30, 31, 32] and references therein for further relevant discussions. 
Such features, and others (e.g. Tully-Fisher relations [33, 34]) suggests to this author that 
dark matter self interactions are likely to be important on small scales, cf. [35, 36, 37]. 

In this paper we consider a particular class of collisional dark matter, dissipative dark 
matter. By this we mean dark matter with strong self interactions which include dissi¬ 
pative particle processes involving the emission of a massless (or very light) dark bosonic 
particle. That is, dark matter is envisaged to have broadly similar particle properties to 
ordinary matter. Multicomponent hidden sector dark matter featuring a massless ‘dark 
photon’ (from an unbroken dark U{1) gauge interaction) which kinetically mixes with the 
ordinary photon provides a concrete example of such dark matter [38, 39, 40, 41] (see 
also the review [42] for a more detailed bibliography). ^ The kinetic mixing interaction 
[52, 53, 38], 


C.nt = , ( 1 ) 

facilitates halo heating by enabling ordinary supernovae to be a source of these ‘dark 
photons’ [54, 55]. These dark photons can transport a large fraction of a supernova’s core 
collapse energy to the halo (potentially up to ~ 10^^ ergs per supernova for kinetic mixing 
of strength e ~ 10“®). Dark matter halos can thereby expand and contract in response 
to these heating and cooling processes, but at the current epoch could (typically) have 
evolved to a steady state or ‘equilibrium’ conhguration where heating and cooling rates 
locally balance. 

The resulting dynamics allows the dark matter density prohle to be related to the 
distribution of ordinary supernovae in the disk of a given galaxy. In a previous paper 
[56] a simple formula was derived encoding this relation which, it was argued, should ap¬ 
proximately represent dissipative dynamics independently of the details of the particular 
dissipative particle physics model. This formula is highly predictive as all of the model 
dependence, fundamental physics etc., was able to be condensed into a single parameter, 
A. It was shown that this approach could potentially explain the apparent correlations 

^More generally, hidden sector dark matter models with an unbroken dark U{1) gauge interaction 
have been quite widely discussed in the literature, e.g. [43, 44, 45, 46, 47, 48, 49, 50, 51]. 


1 



between the dark matter and ordinary matter distributions in galaxies. The purpose of 
the present paper is to further explore this predicted halo prohle of dissipative dark mat¬ 
ter, examining all 26 dwarf galaxies in the LITTLE THINGS sample [19]. In doing so, we 
aim to improve on previous work by modelling the supernovae distribution via the mea¬ 
sured UV and Ha fluxes. This is a more direct measure of the supernovae distribution 
than that considered previously in [56] where a Kennicutt-Schmidt type relation [57, 58] 
was used to model the supernovae distribution in terms of the baryonic gas density. We 
also briefly consider THINGS spirals [59]. Although spirals are generally not as useful in 
testing halo prohles due to the much larger baryonic contribution, they are important, 
especially for exploring the scaling behaviour of A. 

2 Halo profile from dissipative dark matter 

It has been known for some time that dissipative dark matter models can explain dark 
matter phenomena on large scales: LSS and GMB anisotropies [60, 61, 62]. These types 
of models are also consistent with other cosmological probes, such as SNeff{BBN) and 
SNeff{CMB) [63, 41]. The small scale structure of dark matter, especially the dark 
matter distribution around galaxies, has been pursued in [39, 64, 65, 66, 67, 41, 42, 56]. 

In this picture the halo of disk galaxies is in the form of a strongly self interacting gas 
of particles. In the model of [41] this gas is a plasma composed of dark electrons {Fi) 
and dark ions {F 2 ) interacting via massless dark photons, 'jn. Such a halo is dynamical 
and can be modelled as a fluid governed by Euler’s equations, with both heating and 
cooling processes. As in the previous studies, we assume the existence of a kinetic mixing 
interaction, Eq.(l), of strength e ~ 10“®. This small interaction results in a substantial 
halo heat source, with up to around half of the total core-collapse energy of ordinary 
supernovae converted into dark photons [54, 55, 64]. These dark photons propagate out 
into the halo where they can eventually be absorbed there via some interaction process, 
with dark photoionization a likely suspect in the specihc models studied [39, 64, 65, 66, 41]. 
It follows that the heating rate at a particular point, P, in the halo is proportional to the 
product of the dark matter density and dark photon energy flux at that point: r/jeat(r) oc 
n{r)F^jy{r). The halo is dissipative, and cools via dark bremsstrahlung (and potentially 
other processes), which means that the cooling rate at the point, P, is proportional to the 
square of dark matter density: TcooKi') cc n(r)^. [The proportionality coefficients depend 
on the details of the dissipative particle physics, and do not need to be specihed for the 
purposes of model independent considerations.] 

Given sufficient time, a dark matter halo can evolve to a steady-state conhguration 
where the heating and cooling rates balance at every location, Theati'^) = TcooK^)) so that: 

n{r) cc F^^{r). ( 2 ) 

The timescale of this halo evolution, r, is presumed to be much less than the current 
age of the Universe. A rough estimate for this timescale is the time for which dissipative 
interactions would dissipate all the halo’s energy in the absence of heating, i.e. r ~ 
{3/2)n{r)T/Tcooi (here T is the halo temperature). The precise value of r is of course 
model dependent, but expected to be less than around a Gyr in the specihc models 
studied [41, 42]. 

As discussed earlier, the hux of dark photons which govern this dynamics is expected 
to originate within and around core-collapse supernovae. Supernovae are the hnal evo¬ 
lutionary stages of large (M ~ Sm©) stars which are located in the galactic disk. Using 
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spherical coordinates and setting the presumed thin disk at 0 = 7r/2 (i.e. the axis is nor¬ 
mal to the plane of the disk) the energy flux of these dark photons at a point P = (r, 6, 0) 
within an optically thin halo is given by: ^ 


oc 


dcj) dr r 


^SN 




47r[r^ + r^ — 2rr sin 6 cos(0 — 


( 3 ) 


Here, T,sn(j, 0) is the type II supernova rate per unit area in the disk. Strictly, T^sn is 
some appropriately weighted average over the timescale ~ r, expected to be many orders 
of magnitude greater than the time between the discrete supernova events. 

The dynamically driven balancing of heating and cooling rates [r/ieat(r) = TcooK'')] 
dictates the dark matter density, via Eq.(2). That is. 


p(r, 9, (j)) = \ dcj) dr r 


^SN 




47r[r^ +r^ — 2rr sin 6 cos(0 — 0)] 


( 4 ) 


The proportionality coefficient. A, depends on the dark photoionization cross section, su¬ 
pernovae dark photon energy spectrum, and halo properties: ionization state, composition 
and temperature. Evidently, there is a signihcant parameter degeneracy which potentially 
makes this kind of dark matter model extremely predictive as far as galaxy dynamics is 
concerned, despite the nontrivial nature of dissipative dark matter. 

A critical assumption in deriving Eq.(4) is that the dark halo is able to dynamically 
evolve until it reaches a steady-state conhguration. Naturally, this could only occur if the 
galaxy is sufficiently isolated and the heat source (supernova rate) is sufficiently stable. 
Environmental effects, such as mergers, perturbations from nearby galaxies etc. could 
hinder an actual galaxy from attaining this ‘equilibrium’ conhguration. The degree to 
which a particular galaxy is perturbed away from equilibrium would depend on the size 
of the perturbation as well as when it occurred. Once a perturbation ceases, the galaxy 
can re-equilibrate on the timescale ~ r, which as mentioned above is expected to be less 
than a Gyr. An important class of ‘perturbed’ galaxies are those currently in a starburst 
phase. Starburst galaxies exhibit a rapidly varying star formation rate over the last few 
hundred million years. Given that ordinary supernovae are the presumed source of the 
dark photons which heat the halo, this then implies a rapidly changing halo heating rate 
[r/ieat(r, t)]- Such a rapidly changing halo heat source can, of course, potentially perturb 
the halo away from the steady-state configuration. Thus, the density prohle of the dark 
matter halo of starburst galaxies can depart signihcantly from Eq.(4). 

Although this discussion has been couched in terms of a specific kind of dissipative 
dark matter model (where halo heating is supplied by massless dark photons produced via 
kinetic mixing induced processes in ordinary supernovae) Eq. (4) can hold more generally. 
The basic assumptions are only that dark matter is dissipative, behaves as a fluid which 
evolves to a steady-state conhguration where heating and cooling rates equilibrate (on a 
timescale r ~ 1 Gyr). It is easy therefore to envisage many other dissipative dark matter 
models leading to Eq.(4). For example, models where heating is transported to the halo 

^ Supernovae produce a flux of dark photons with an uncertain energy spectrum. In general the optical 
depth is a frequency (energy) dependent quantity so that in actuality the halo might only be optically 
thin for a range of dark photon frequencies. However these optically thin dark photons can dominate 
the energy transport, making the optically thin approximation potentially valid even in this case. This 
conclusion is also supported by the numerical work of [65, 66] where the effects of finite optical depth 
was considered where the uncertain dark photon supernova spectrum was modelled with a wide range of 
possible dark photon frequency distributions. 
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via massive dark photons are possible. Alternatively heating might be transported from 
ordinary supernovae to the halo via light scalar particles. In these scenarios, various 
considerations constrain the mass of the massive dark photon or dark scalar to be around 
~ 10 keV. However, models with energy sources other than supernovae, e.g. ordinary main 
sequence stars, are unlikely given the typically stringent constraints on exotic energy loss 
mechanisms for such objects. 

In general, the coefficient A in Eq.(4) will depend on the position, r, 6*, 0 and can 
be model dependent. Such spatial variation of A is due, in part, to the dependence on 
A on the halo temperature. Simplihed calculations, though, within a given dissipative 
model [65, 42] suggest that the spatial variation of the temperature is only important 
in the central region of the galaxy and even in that region can (typically) be relatively 
modest for small galaxies. For the purposes of the analysis given here we shall consider 
A in Eq.(4) to be spatially independent as a zeroth order approximation. Under this 
approximation the spatial dependence of p(r, 9, 0), Eq.(4), is determined solely by that of 
the flux, F^jy{r,9,(l)). 

Eq.(4) is subject to another signihcant caveat. It provides the mass density only of the 
diffuse dark matter (plasma) component. Dissipative dark matter models can also have 
dark matter in the form of ‘dark stars’, that is, stars composed of dark matter particles. 
Naturally it is quite challenging to figure out the proportion and distribution of such a 
component, and for the present discussion it will be assumed negligible. ^ 

With the above important qualihcations, Eq.(4) provides a description of the dark 
matter density in a given galaxy. Disk galaxies though come in various sizes, how does 
A depend on the galaxy scale? The physical quantity of most interest is again the halo 
temperature, T. Assuming that the halo is isothermal and in approximate hydrostatic 
equilibrium, the temperature can be expressed in terms of the maximum value of the 
rotational velocity: [39, 41, 42] 


T = \m (5) 

where rh = '^Uimi/ '^Ui is the mean mass of the particles making up the dark plasma. 
The maximum value of the rotational velocity, varies between around 20 km/s for the 
smallest dwarfs to around 300 km/s for the largest spirals. If we assume that the heating 
is due to photoionization of K-shell dark atomic states and that the binding energy of 
these states is larger than the halo temperature of the largest spirals, then the heating 
rate is approximately independent of the halo’s temperature. The dark photons heating 
the halo encounter fully occupied K-shell states for all galaxy scales of interest. Cooling, 
on the other hand, is expected to depend more strongly on the halo’s temperature since 
the bremsstrahlung cross section depends on temperature, leading to: Tcooi oc \/Tn(r)^. 
Thus, these simplihed dissipative models suggest that A oc 1/\/T, and via Eq.(5) give 
A oc l/n™"’”. However more complicated behaviour is certainly possible, and would be 
expected for temperature regions where other cooling processes become important, ^ the 

^Note though that for elliptical galaxies (and possibly dwarf spheroidals) the situation is expected to 
be very different. For these galaxy types there is (currently) very low star formation rate and therefore 
heating of the dark halo from ordinary supernovae is not expected to be significant. Any diffuse dissipative 
dark matter halo around these kinds of galaxies would therefore have collapsed into ‘dark stars’ [42, 41]. 

®For an important class of dissipative models, including those studied in [39, 64, 65, 66, 42, 41], dark 
atom recombination and line emission can potentially significantly contribute to cooling. These processes 
also feature Tcooi oc n(r)^ and thus Eq.(4) still holds, but the proportionality constant can have a more 
complicated dependence on halo temperature. If such processes are dynamically important they can 
modify the predicted \{T) oc Xj'jT scaling. 
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details will of course depend on the particle physics properties of the specihc dissipative 
model. As will be discussed in more detail in section 3, this galaxy scaling relation is 
related to the Tully-Fisher relation [33]. 

The halo dark matter contribution to the gravitational acceleration at a point in the 
plane of the disk can be straightforwardly calculated given the density prohle, Eq.(4), and 
Newton’s law of gravity. Assnming an azimnthally symmetric disk, the motion of the gas 
and stars in the disk is circular with speed Vhaio given by: 


Vhalo 

r 


Gn dcj) dcos6 dr r 


.2 p{r, 6) cos a; 




( 6 ) 


Here, d? = — 2rf sm.6 coscj)^ cos a; = (r — r sin^cos^)/^ and is Newton’s 

constant. 

In the previous paper [56] the density prohle, Eq.(4), was examined by hrst modelling 
the snpernovae rate with an exponential disk: T^sNir) = (i?5Ar/27rr|,) [68]. This 

could only be a very crude approximation, since the exponential disk is an approximate 
measnre of the stellar population, including middle aged and older stars, while the qnantity 
we reqnire is the recent (large) star formation rate. Anyway, if S 5 Ar(r) is an exponential 
fnnction of radius then Eq.(4) becomes ronghly eqnivalent (as far as rotation curves are 
concerned) to the cored distribntion: pisoir) = Po'f'o/i'^^ + with tq ^ ro- Snch 
a constrained quasi-isothermal prohle is known to be phenomenologically snccessfnl in 
explaining the observed shapes of rotation cnrves, e.g. [15, 69, 70, 71, 30, 19]. Dissipative 
dark matter therefore provides an underlying theoretical explanation for the snccessfnl 
qnasi-isothermal prohle. 

The dark matter density prohle of Eq.(4) is, in fact, further constrained. At small 
radii, r < r^), the dark matter density [Eq.(4)] can be related to the central surface 
density, S 5 Ar( 0 ) = RsN/{‘27Trjj), via: 


p(r) 


ASgAr(O) 

2 



+ constant 


(7) 


This implies a scaling relation connecting the inner circnlar velocity gradient with this 
central surface density: v oc r^J \T,sn{^)- To connect with measnrable qnantities, one 
could replace Ss'7v(0) with the central snrface brightness in the UV band, S;7y(0) and 
neglect the variation of A (which is anticipated to be fairly weak for dwarf irregnlar 
galaxies). This leads to the rongh scaling relation: v oc r^T^uviS^) or 


log 



—0.2/io + constant 


( 8 ) 


where /io is the central surface brightness in magnitnde nnits. Such a scaling relation has 
in fact been observed to hold for spiral and irregnlar galaxies [32]. 

In addition to modelling T^sn via an exponential disk, the previous paper [56] further 
considered modelling in terms of the baryonic gas density via a Kennicntt-Schmidt 
type relation [57, 58]. This allowed a connection between the supernovae distribntion 
with current local properties of a given galaxy. In the following we shall consider a more 
direct estimate of the snpernovae distribntion, by modelling the snpernovae distribntion 
via the measnred UV and Ha flnxes. 
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3 Dissipative dark matter versus dwarf galaxies 
3.1 The LITTLE THINGS dwarfs 


In this paper we shall examine all 26 dwarf galaxies comprising the LITTLE THINGS 
sample [19]. An inspection of these 26 galaxies suggests a loose classihcation: 18 dwarfs 
which feature ‘classically’ shaped rotation curves with a (typically) linear rise in Vrot near 
r = 0 smoothly transiting to a flat rotation curve at the greatest measured radii, three 
dwarfs which show a classically shaped rotation curve for r less than some radius, R*, but 
show a downturn for r > R*, and four dwarfs which show a classically shaped rotation 
curve for r less than some radius, i?i, but feature a ‘hump’ at r > i?i. In addition there 
is one dwarf (DD047) whose rotation curve is irregularly shaped. In figure 1 we give an 
example for each of these galaxy types. 

The theoretically predicted rotation velocity, is the sum of the various contribu¬ 
tions added in quadrature: 


Vrot = + V^as + V%ars ■ (9) 

Here, Vgas and Vstars are the baryonic gas and stellar contributions while Vhaio is the dark 
matter contribution, Eq.(6). The latter is given in terms of the density, Eq.(4), which 
can be evaluated if the supernovae distribution of the galaxy is known. To proceed we 
therefore need to model the supernovae rate, T,sn{v). 
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Figure 1: Examples of the four classes of galaxies in the LITTLE THINGS sample; DD052, 
classically shaped rotation curve (linear rise, transiting to flat), F564-V3, classically shaped but 
with a downturn in Vrot, DDO70, classically shaped but with a ‘hump’ at large R, DD047, an 
irregularly shaped rotation curve. 
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In the literature UV and Ha emission have both been used as traces for star formation 
rates, see e.g. [72, 73, 74], The UV radiation is directly emitted from the photospheres 
of large stars, with M* ~ 3Mq (O- through later-type B-stars) while the Ha nebular 
emission results from the recombination of hydrogen which can only be ionized by very 
large stars, with M* ~ 17M0 (the most massive O- and early- type B-stars). The Ha flux 
provides a measure of star formation rate on fairly short timescales of ~ 10 Myr, while 
UV flux provides a measure on a somewhat larger timescales: ~ 100 Myr. The timescale 
most relevant for the supernovae rate, T^sn, in Eq-(4) is some appropriately weighted av¬ 
erage over the timescale characteristic of the dissipative dynamics, r ~ {3/2)n{r)T/Tcooi, 
expected to be less than around a Gyr as mentioned earlier. A priori, both UV and Ha 
fluxes present reasonable candidates for the relevant supernovae rate, and we will consider 
both prohles in the analysis to follow. 

The Galaxy Evolution Explorer (GALEX) [75] has obtained images of nearby galaxies 
in the far UV (FUV) bandpass of 1350-1750 A and near UV (NUV) bandpass of 1750-2800 
A. The GALEX data was analyzed in [76] for a large selection of dwarf galaxies, including 
24 of the LITTLE THINGS sample of 26 galaxies, yielding the azimuthally averaged 
surface brightnesses, HfuvA)- This quantity can be converted to the flux density via 
the standard formula: FfuvA) — erg cm“^ s“^ Hz“^ arcsec”^, and 

similarly for the NUV flux. Also given in that reference is the Ha surface brightnesses, 
Ha, which can also be converted to a Fuaij') flux. In the following, we shall model 'LsnA) 
with each of these fluxes. Of most interest is the radial dependence of these quantities as 
the overall proportionality constant can be absorbed into a redefinition of the coefficient 
A of Eq.(4); this redehnition is denoted as A (and will be dehned more precisely later on). 

Gonsider Erst the 18 classical dwarfs and DD047. In table 1 we give some properties 
of these galaxies: the distance of the galaxy as given in [19] and also the absolute AB 
FUV-band magnitude (evaluated from the apparent magnitude assuming this distance 
measurement). ® Also given in table 1 is = x^/Ffdof (for a 1 parameter fit, the number 
of degrees of freedom is N^of = N^ata — 1) for the fit of the model to the rotation curve 
data, modelling T^snA) with a) FfuvA')) b) F^uvA) and c) Fua, each obtained from [76] 
as described above. ^ Here, the yA function for each galaxy is constructed in the usual 
way: 

^data 

E 

i=l 

where the N^ata (asymmetric drift corrected) binned rotation curve measurements, (i)± 
Svlotii), are obtained from [19] and Vrot is evaluated via Eq.(9) (with Vgas and vstars from 
[19]). This function was minimized with respect to variations of the parameter A 
(which was allowed to have a different value for each galaxy). The related quantity. A, to 
be defined in Eq.(14), is also given in the table. 

The values given in the table suggests that the density profile, Eq.(4), motivated 
by the dissipative dynamics is consistent with a zeroth order approximation for most of 
the 18 classical dwarfs. The resulting fits are shown in figures 2-4 for the 18 ‘classical’ 
dwarfs and DD047. In these figures, the FUV flux was used as the tracer for 

®For dwarf galaxies extinction corrections are generally small. For the estimated FUV absolute magni¬ 
tude of the dwarf galaxies the internal and galactic (foreground) extinction corrections were nevertheless 
included following the procedure of [74, 73] using data obtained by the online updated nearby galaxy 
catalogue [77]. 

^For the dwarfs: DD052, DDOlOl, DDO210, DD0216. there were insufficient Ha flux measurements 
for which to construct the flux, while for ICIO UV flux measurements were unavailable. 


'^rot ) 


exp / -x 

VrotA) 


SVrotiA 


( 10 ) 
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very similar results occur for the other two traces considered. Indeed, as the Xr values 
in the table indicate, use of the NUV flux generally leads to a slightly better £t to the 
rotation curve data, while with the Ha flux the fit fares slightly worse. It is also noted 
that use of UV flux generally leads to a slightly better fit to the rotation curves than the 
Kennicutt-Schmidt-type power law considered in [56].® 

The density profile, Eq.(4), is highly constrained as it depends on just one parameter 
A. This parameter sets the normalization of the dark halo contribution of the rotation 
curve. The shape of the curve is fully predicted by Eq. (4) once the supernovae distribution 
is input via the FUV flux radial profile. Figures 2-4 indicate that the radius where the 
rotation curve transitions from linear to flat is generally reproduced reasonably well. ® 
Note though that the rotation curve error bars are not purely statistical as they are 
usually derived by taking the difference between rotation curve measurements from each 
side of the galaxy. It is possible that the small but perceptible difference between the fit 
and the data of some of these galaxies is an indication for corrections to the zeroth order 
approximation of Eq.(4). There are many possible sources of corrections including: radial 
dependence of the halo temperature, modelling of the supernova distribution, correcting 
for azimuthal asymmetry etc. A simple way to incorporate first order corrections is to 
allow A in Eq.(4) to be spatially dependent: A(r) = Aq + Air -|- ... Keeping only the first 

two terms gives a dark matter density defined in terms of two parameters, Aq, Ai. The 

two parameter fit for two examples, DD0154 and WLM is shown in figure 5. As this 
figure illustrates, the inclusion of first order corrections in this simple phenomenological 
manner allows for an improved representation of the data. 

We now consider the three dwarfs with a downturn in Vrot- A possible explanation for 
the downturn is that it represents the physical extent of the dark halo for these galaxies. 
If this is the case then we could try to model this boundary effect with a sharp cut-off: 

A(r) = A for r < /?*, 

A(r) = 0 for r > i?* . (11) 

This introduces one additional parameter, R*. The values for this two parameter 
(A, R*) £t are given in table 2 and the resulting fits are shown in figure 6. It is of course 
also possible that the downturn is an indication that the halos of these dwarfs are not 
currently in a steady-state configuration due to some perturbation. This is also a possible 
explanation of the dwarfs with a ‘hump’, as we will now examine. 

The four remaining dwarfs of the LITTLE THINGS sample are those which feature a 
‘hump’ at the boundary region. The hump feature might be taken as an indication that 
the halos of these dwarfs are not entirely in a steady-state configuration where heating 
and cooling rates locally balance. This situation could conceivably arise if a galaxy was 
being perturbed, either externally by a nearby galaxy or via some internal mechanism. 
Such perturbations, if significant, would affect the star formation rate of the galaxy. This 
suggests that a galaxy’s recent star formation history could possibly be used as a diagnostic 
aid in determining whether its dark halo is likely to be in a steady-state configuration. 

It may also be possible for a galaxy to be in a partially equilibrated state. Perhaps 
only the inner part of the halo could be in an approximate steady-state configuration. 

curious exception is DD047, for which the Kennicutt-Schmidt type power law {Y^snR) oc 
[E gos(r)]'^) leads to a significantly better fit: = 5.2 (for exponent N = 2), ^^{KS) = 6.7 

(for exponent N = 3) cf. Xr(FUV) = 9.7. 

®In fact, this reproduction is also insensitive to the uncertainty in the distance of the galaxy, D. This 
can most easily be seen by working with the angular variable, r/Z?, and noting that both rotation curves 
and surface photometry are measured in terms of such a quantity. 



Galaxy D (Mpc) Mpuv Xr (FUV) Xr (NUV) 

Xr {Ha) A (km^/s^) 

18 classical dwarfs 



DD043 

7.8 

-13.14 

0.55 

0.69 

0.72 

5.32E-3 

DDO50 

3.4 

-15.42 

3.53 

3.52 

3.06 

4.36E-4 

DD052 

10.3 

-13.36 

0.18 

0.15 

- 

1.50E-2 

DD053 

3.6 

-12.51 

1.30 

1.25 

1.30 

4.26E-3 

DD087 

7.7 

-13.09 

0.96 

0.79 

2.54 

1.41E-2 

DDOlOl 

6.4 

-11.59 

28.6 

22.9 

- 

0.11 

DD0126 

4.9 

-13.39 

0.56 

0.55 

0.92 

4.62E-3 

DD0133 

3.5 

-13.00 

5.1 

4.1 

6.4 

1.22E-2 

DD0154 

3.7 

-13.10 

0.81 

0.83 

0.80 

l.OlE-2 

DDO210 

0.9 

-8.23 

0.84 

0.69 

- 

4.53E-2 

DD0216 

1.1 

-9.48 

0.50 

0.53 

- 

2.98E-2 

ICIO 

0.7 

- 

- 

- 

0.07 

- 

NGC1569 

3.4 

-16.80 

0.61 

0.61 

0.50 

9.81E-5 

NGG2366 

3.4 

-15.32 

0.28 

0.27 

0.23 

1.90E-3 

NGG3738 

4.9 

-14.72 

6.14 

4.27 

3.32 

1.56E-2 

WLM 

1.0 

-12.51 

0.53 

0.44 

0.92 

9.29E-3 

Haro29 

5.9 

-13.68 

0.40 

0.40 

0.55 

3.71E-3 

Haro36 

9.3 

-14.61 

1.86 

1.63 

1.90 

1.72E-3 

A dwarf with irregularly shaped rotation curve 

DD047 

5.2 

-13.83 

9.7 

10.8 

12.5 

9.33E-3 


Table 1: LITTE THINGS classical dwarfs (defined in text) and DD047. D {Mfjjv) is the 
distance (FUV absolute AB magnitude) of the dwarf and Xr given for the 1-parameter (A) fit 
of the rotation curve data to the model, Eq.(6), Eq.(4), modelling T^sn with a) EUV, b) NUV 
and c) Ha surface brightness radial profiles. The fitted parameter A, Eq.(14), is also given. 


Galaxy 

D (Mpc) 

Mpuv 

xl with R* (kpc) 

A (kmVs2) 

DD046 

6.1 

-13.15 

3.85 with R* = 1.90 

3.19E-2 

DD0168 

4.3 

-13.64 

3.20 with R* = 3.07 

9.10E-3 

F564-V3 

8.7 

-10.94 

0.56 with R* = 2.15 

5.47E-2 


Table 2: LITTLE THINGS dwarfs with a downturn in the rotational velocity. D (Mpiiv) is 
the distance (FUV absolute AB magnitude) of the dwarf and Xr is given for the 2-parameter 
(A, R*) fit of the rotation curve data to the model, Eq.(6), Eq.(4), Eq.(ll). The supernovae 
distribution, T^sn, is here modelled with the FUV surface brightness radial profile. The fitted 
parameter A, Eq.(14), is also given. 
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Galaxy 

D (Mpc) 

Mpuv 

Xl (FUV) for r <Ri (kpc) 

A (kmVs2) 

GVnldwA 

3.6 

-11.56 

1.24 for r < 1.4 

1.41E-3 

DDO70 

1.3 

-11.92 

0.71 for r < 1.6 

1.19E-2 

IG1613 

0.7 

-12.80 

1.15 for 0.3 < r < 1.7 

2.38E-4 

UGG8508 

2.6 

- 

1.72 for r < 1.4 

- 


Table 3: LITTLE THINGS dwarfs with a ‘hump’ in the rotational velocity. D (Mpjjy) is 
the distance (FUV absolute AB magnitude) and Xr is given for the 1-parameter (A) fit of the 
rotation curve data to the model, Eq.(6), Eq.(4), for the inner radial region, r < Ri. The 
supernovae distribution, Tisn-, is here modelled with the FUV surface brightness radial profile, 
except for UGC8508 which was modelled with Ha (FUV was unavailable in the case). The 
fitted parameter A, Eq.(14), is also given. 


This situation might arise because of the disparate timescales: the higher density inner 
region of a dark halo can cool and equilibrate more rapidly than the less dense outer 
region. The dwarfs with a ‘hump’ (and possibly also those with a ‘bump’) might result 
from such a condition. In table 3 and hgure 7 we show the results of the £t to the inner 
region of the dwarfs with a ‘hump’. 

It is perhaps noteworthy that some of the dwarfs which feature poor fits to the steady- 
state solution, Eq.(4), are known to be undergoing starburst activity. In particular, of the 
four dwarfs with the lowest A values, DDO50, NGC1569, IC1613 and CVnldwA, three 
of them [NGC1569, DDO50 (Holmberg II) and GVnldwA (UGGA292)] are known to be 
in a starburst phase, e.g. [74, 78]. In fact, a low value of A is an indication that halo 
heating exceeds cooling {Theat > ^cooi) a situation which would imply that these halos are 
currently expanding. 

Dissipative dynamical halos of the type envisaged here would be expected to play an 
important role in regulating the star formation rate of a given galaxy. Any perturbation 
either external or internal that causes the halo to move out of equilibrium, Theat r cool 1 
would lead to an expanding or contracting halo dynamically adjusting under the various 
forces. Such dynamics, which could in principle be described by Euler’s equations of 
fluid dynamics, would be expected to be a major influence on ordinary star formation 
rates. As the halo expands or contracts, the baryonic gas density decreases or increases in 
response to the changing gravity. These physical changes in baryonic gas density would 
in turn correlate with the star formation rate (as suggested by the Kennicutt-Schmidt 
relation). This interplay between the dynamical halo and star formation rates could be 
one of the main drivers of starburst activity. It is a temporary non-equilibrium phase 
where Tcooi 7^ Theat and the halo is undergoing significant changes influencing the star 
formation rate and indeed vice versa. 

The dwarf DDOlOl has the distinction of featuring the worst fit, with Xr{FUV) = 
28.6. It was emphasised in [19] that the inner slope of DDOlOl does have potentially large 
uncertainty due to possible beam smearing effects. Whether such effects could reconcile 
the predicted rotation curve with the data is unclear. If this (or other) systematic effects 
are not responsible for the disagreement, then one could speculate that this very faint 
dwarf has evolved to the point where there is insufficient heating to support the halo. 
The halo may have cooled and collapsed (or is in the process of collapsing) into dark 
stars. Indeed table 1 indicates that the A value for this dwarf is anomalously large, which 
is evidence in favour of this interpretation. 
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Figure 2: Rotation curves for the dwarf galaxies, DD043, DDO50, DD052, DD053, DD087, 
DDOlOl, DD0126 and DD0133. The data is from [19]. The stellar (dashed line), baryonic gas 
(dotted line) and dissipative dark halo (dashed dotted line) contributions are shown. The solid 
line is the sum of these contributions added in quadrature. The halo contribution results from 
the 1-parameter fit, Eq.(4) with A assumed spatially independent, [^sn in Eq.(4) is modelled 
via the FUV surface brightness radial profile from [76] as described in the text.] 
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Figure 4; Rotation curves for the remaining two ‘classically’ shaped dwarfs: Haro29 and Haro 
36 and also the irregularly shaped rotation curve of DD047. Notation as in figure 2. 




Figure 5; Rotation curves for DD0154 and WLM. Notation as in figure 2, except that the halo 
contribution was calculated from Eq.(4) allowing for a spatially dependent A; A(r) = Aq + Air 
and thus a two-parameter fit. 
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Figure 6: Rotation curves for the three dwarf galaxies, DD046, DD0168, F564-V3, indicating 
a downturn in Vrot- The downturn is modelled with a hardcut off in the dark matter density, as 
described in the text. Notation as in figure 2. 






Figure 7; Rotation curves for the remaining LITTLE THINGS dwarfs [CVnldwA, DDO70, 
IC1613, UGC8508] which feature a ‘hump’ at the largest measured radii. Shown are fits to the 
inner part of the rotation curve data, see discussion in text. Notation as in figure 2. 
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Galaxy D (Mpc) Mpuv A (km^/s^) 


NGC925 

9.2 

NGC2403 

3.2 

NGC2841 

14.1 

NGC2903 

8.9 

NGC2976 

3.6 

NGC3031 

3.6 

NGC3198 

13.8 

NGC3521 

10.7 

NGC3621 

6.6 

NGC4736 

4.7 

NGC5055 

10.1 

NGC6946 

5.9 

NGC7331 

14.7 

NGC7793 

3.9 


-18.28 

6.1E-4 

-17.85 

1.2E-3 

-18.77 

2.0E-3 

-18.45 

1.2E-3 

-14.67 

6.6E-3 

-18.08 

4.8E-4 

-18.97 

4.4E-4 

-18.37 

1.5E-3 

-18.29 

8.8E-4 

-16.86 

6.5E-4 

-18.57 

9.4E-4 

-18.74 

8.4E-4 

-18.96 

7.5E-4 

-17.18 

1.5E-3 


Table 4: THINGS spirals; D is the distance of the spiral and Mpjjv is the FUV absolute 
AB magnitude corrected for internal and foreground extinction. Also given is the parameter A, 
Eq.(14). 

3.2 THINGS spirals 

Let us now briefly comment on spiral galaxies. Compared to dwarf irregular galaxies, 
spirals have a much larger stellar component, which typically dominates over the dark 
matter contribution in the inner part of the galaxy. This stellar component is poorly 
constrained due to uncertain stellar mass-to-light ratio, T*. This unfortunately tends 
to make spirals much less useful in testing dark matter models. Also, the UV and Ha 
flux measurements for spirals have larger uncertainties due to much larger extinction 
corrections. A sample of spirals with high resolution rotation curves has been given 
by the THINGS collaboration [59]. An example is given in hgure 8, where the £t to 
the rotation curve of NGC2403 is given. In this example the baryonic contribution was 
obtained from [59] and modelled allowing T* to vary, and the dark matter contribution 
was obtained from Eq.(6) and Eq.(4). The supernovae rate in Eq.(4) was modelled via a 
Kennicutt-Schmidt relation, with exponent N = 2, T^sn oc (UV surface profile was 
unavailable). 
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Figure 8: The spiral galaxy, NGC2403. The stellar (dashed line), baryonic gas (dotted line) and 
dissipative dark halo (dashed dotted line) contributions are shown. The solid line is the sum of 
these contributions added in quadrature. 


3.3 The scaling of A 

In the fits of the dwarf galaxies the parameter A in Eq.(4) was allowed to have a different 
value for each galaxy. In specific dissipative dark matter models this parameter is expected 
to have a more universal character, with simple arguments suggesting a weak scaling: 
A oc [see earlier discussion around Eq.(5)]. 

At sufficiently large radii the dark photon flux falls at the geometric rate: oc 1/r^. 

At such large radii the dark matter halo mass density profile [Eq.(4)] reduces to: 


p(r) 


. Rsn 


( 12 ) 


Here Rsn is the integrated supernova rate in the disk. Eq.(6) gives an asymptotically flat 
value for the halo contribution to the circular velocity: 


•2 = \GnR 


^halo 


NJ^SN 


(13) 


Supernovae are the final evolutionary stage of large stars. As discussed earlier, the FUV 
luminosity {Lfuv) is expected to provide a reasonably direct estimate of the star formation 
rate, so that: Rsn oc Lpuv oc ig convenient then to introduce the quantity 

A: 


A = 


V 


2 

halo 


\Q-0AMfuy 


(14) 


With this definition Eq.(13) implies A= XGnCi where ci = For each 

galaxy, Vhaio can be obtained from the £t of the model to the rotation curve data (for 
the dwarfs with a downturn in Vrot-, listed in table 2, we set Vhaio = The absolute 

FUV AB magnitude, Mpjjvi has already been provided in tables 1-3 for the 24 dwarfs 
for which FUV flux measurements are available. These tables also give the resulting A 
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values, as defined in Eq.(14), for each dwarf. We have also estimated the A values for the 
THINGS spirals. For these spirals we obtained Vhaio from the quasi-isothermal £t with 
free mass-to-light stellar mass parameter as given in [59]. The absolute FUV magnitude 
was obtained using data from [77, 79] and corrected for both internal and foreground 
extinction following the procedure of [74, 73]. The results of this exercise are given in 
table 4. Figure 9 summarizes these results. 

Under the assumptions that the halo heating is approximately independent of the 
halo temperature, T, and halo cooling is dominated by bremsstrahlung, Tcooi oc \/Tn^, 
one expects the rough scaling A oc 1/\/T. As discussed earlier in section 2 [see discussion 
around Fq.(5)], this leads to the correlation: A oc Assuming that ci (defined 

above) is constant this then implies A oc or equivalently, 

Lfuv oc • (15) 

This scaling is also shown in figure 9. For many galaxies one could make the rough 
approximation: Vhaio ~ i from which Fq.(15) reduces to: Lpuv oc [n™“^]^. For spiral 

galaxies, this is approximately equivalent to the Tully-Fisher relation [33], e.g. the study 
[80] found that Lb oc [n)()“*]", where a = 3.4 ±0.09 (the difference between this exponent 
and the value a = 3 is consistent with the slight scaling difference between the FUV and 
B-bands). 

Figure 9 shows a significant scatter of the A parameter for the LITTLF THINGS dwarf 
galaxies. Some scatter can be due to uncertainties, e.g. the determination of the galaxy’s 
absolute magnitude will be affected by the uncertainty in its distance, the parameter c\ 
might vary somewhat between different galaxies etc., however such uncertainties are un¬ 
likely to be the whole explanation The size of the scatter might indicate fundamental 
physical differences between galaxies, or it could simply be an indication that the dissipa¬ 
tive halos are not all in a steady-state configuration. Halo heating and cooling rates may 
be unbalanced due to perturbations from nearby galaxies, current (or recent) starburst 
activity etc. As already discussed, such activity might cause distortions including, possi¬ 
bly, the ‘bump’ and ‘hump’ in the rotation curves apparent in a subset of the dwarfs. If so, 
then it could be a factor behind the poor fits to the shape of the rotation curve that some 
of the dwarfs feature. In figure 10 we therefore consider only a subset of galaxies: the 
THINGS spirals together with the LITTLE THINGS classical dwarfs and only those with 
XriFUV) < 2. This means that the dwarfs DDOlOl, DDO50, DD0133, NGG3738 (and 
those of tables 2 and 3) are excluded, and we also exclude NGG1569 which is undergoing 
extreme starburst activity [78] (whose halo is therefore unlikely to be in an equilibrium 
configuration). 

Figure 10 shows a much reduced scatter (cf. figure 9) of the A values with respect 
to the A oc l/u™“‘^ scaling. This is some evidence in favour of a ‘universal’ A parameter, 
presumably set in part by the fundamental particle physics of the dissipative model. 
Further studies, with a larger sample would be a useful next step. In particular, the 
possibile validity of Eq.(15) over the entire range from dwarfs to spirals warrants further 
investigation, along with possible connection with the baryonic Tully-Fisher relation [34]. 

^*^The uncertainty in the parameter A was estimated from the distance uncertainty, Vhaio uncertainty, 
and the uncertainty in the extinction corrections. The distance uncertainty for most of the galaxies is 
small, typically 10%, while the Vhaio uncertainty was taken from the rotation velocity data points at large 
radii. The uncertainty in extinction was assumed to be ±20% of the extinction correction. 

^^Note however that since Eq.(4) defines the dark matter density directly in terms of the supernova 

density (rather than say, baryonic mass), the estimated values for A should be independent of the uncertain 
mass-to-light ratio of the galaxy. 
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Figure 9: The parameter A [Eq.(14)] versus for the LITTLE THINGS dwarfs (squares and 
circles) and THINGS spirals (stars). The squares are the classical dwarfs while the circles are 
all the others. The solid line is the A oc scaling discussed in text. 



Figure 10: Same as figure 9 except only THINGS spirals and ‘best fit’ classical dwarfs are 
included (see text). 
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4 Conclusion 


There is substantial evidence from rotation curves that dark matter halos have nontrivial 
dynamics. Most striking is that the observed rotation curves are strongly correlated with 
the baryon content. The Tully-Fisher relation provides a global relation, while the dark 
matter rotational velocity generally transits from linear rise (near the galactic center) to 
flat at a distance scale associated with the baryonic distribution, and there are also many 
other correlations. 

Dark matter halos around disk galaxies can have snch nontrivial dynamics if dark 
matter is strongly self interacting and dissipative. Mnlticomponent hidden sector dark 
matter featuring a massless ‘dark photon’ (from an nnbroken dark U{1) gange interaction) 
which kinetically mixes with the ordinary photon provides a concrete example of such 
dark matter. The kinetic mixing interaction facilitates halo heating by enabling ordinary 
snpernovae to be a sonrce of these dark photons. Dark matter halos can expand and 
contract in response to the heating and cooling processes, but for a sufficiently isolated 
and unperturbed halo should have evolved to a steady state or ‘equilibrium’ conhguration 
where heating and cooling rates locally balance. This dynamics allows the cnrrent dark 
matter density prohle to be related to the distribution of ordinary snpernovae in the disk 
of a given galaxy, Eq.(4). Natnrally, this is a simplihed description and would likely be 
invalid for galaxies strongly influenced by pertnrbations, starburst activity etc. 

The ordinary snpernovae distribntion was here modelled via the UV light emission 
prohle, which should provide a fairly direct measurement of this distribution in a given 
galaxy. This improves upon previous work of [56] which used the gas density via a 
Kennicutt-Schmidt type relation. The resulting halo rotation curve then depends on only 
one parameter A (assnmed to be spatially independent as a zeroth order approximation) 
and tested against the full LITTLE THINGS sample of 26 dwarf galaxies. The analysis 
indicates that for most of the LITTLE THINGS dwarfs, the dark matter predictions 
resnlting from the simple formula, Eq.(4), agree reasonably well with the observations. 
Not all of the dwarfs analysed have good agreement, however the poor hts of some of 
these dwarfs can plausibly be ascribed to environmental inhnences which pertnrb the 
halo. Such a perturbation may only be temporary, the timescale depending on the natnre 
of the perturbation and the details of the dissipative particle physics. 

In the simplihed approach considered here, the steady-state conhgnration, Eq.(4), is 
described by a single parameter, A. This parameter is a function of various fundamental 
parameters, e.g. in the model of [41] one has dark particle masses, dark hne strnctnre 
constant etc. The parameter A also depends on less fundamental parameters, including 
galaxy-dependent quantities such as halo temperature. Although in principle A can be a 
complicated fnnction of these parameters, some simple argnments favour an approximate 
scaling with the halo temperature, A oc l/y/T. Such a simple scaling is roughly consistent 
with the analysis of the LITTLE THINGS dwarfs and THINGS spirals presented here. 
Deviations from this scaling are possible and can in principle be calculated within a given 
dissipative particle model, and could thereby provide an important means of testing and 
constraining the underlying particle physics. 
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